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In the collisional region at finite temperatures, the collective modes of superfluids are described by 
the Landau two-fluid hydrodynamic equations. This region can now be probed over the entire BCS- 
BEC crossover in trapped Fermi superfluids with a Feshbach resonance, including the unitarity 
region. Building on the approach initiated by Zaremba, Nikuni, and Griffin in 1999 for trapped 
atomic Bose gases, we present a new variational formulation of two-fluid hydrodynamic collective 
modes based on the work of Zilsel in 1950 developed for superfluid helium. Assuming a simple 
variational ansatz for the superfluid and normal fluid velocities, the frequencies of the hydrodynamic 
modes are given by solutions of coupled algebraic equations, with constants only involving spatial 
integrals over various equilibrium thermodynamic derivatives. This variational approach is both 
simpler and more physical than a direct attempt to solve the Landau two-fluid differential equations. 
Our two-fluid results are shown to reduce to those of Pitaevskii and Stringari for a pure superfluid 
at T = 0. 

PACS numbers: 03.75.Kk, 03.75.Ss, 67.40.-w 



I. INTRODUCTION 

In 1938, Tisza pj first suggested that the classic fea- 
tures of superfluid helium were a manifestation of a two- 
fluid hydrodynamics, originating from a Bose-Einstein 
condensate. The correct two-fluid equations in the non- 
dissipative limit were formulated and solved by Landau a 
few years later [U-Q. Two- fluid hydrodynamics describes 
the coupled dynamics of the superfluid and normal fluid 
components when collisions produce a state of local ther- 
modynamic equilibrium. This collision-induced state is 
described by the usual thermodynamic variables, which 
are now dependent on position and time. The two-fluid 
equations are essentially the same as the fluid dynamics 
of a normal fluid, but with additional equations describ- 
ing the new superfluid component. In simple terms they 
describe the collective oscillations of frequency to such 
that lvtr <c 1, where tr describes the time it takes for a 
nonequilibrium state to reach local equilibrium (see, for 
example, Ref. 0). For brief accounts of two- fluid hydro- 
dynamics in the context of uniform quantum gases, see 
Refs. H||. 

The two-fluid regime was difficult to reach in the first 
wave of experiments on trapped Bose-condensed gases 
after 1995. The typical collisional cross-section and the 
achievable densities were not sufficient to reach the region 
described by Landau's two-fluid hydrodynamics. An ex- 
ception was the pioneering experiment by Ketterle and 
coworkers 0, which later theoretical work showed was 
well within the two- fluid region 4]. However, there is 
new interest in the study of the two-fluid collisional dy- 
namics in trapped gases. One reason is that with atom 
chips, one can now produce tightly localized high den- 
sity atomic Bose condensates. In addition, in a trapped 
two-component Fermi gas, one can use a Feshbach res- 
onance to adjust the magnitude and sign of the s-wave 
scattering length a between Fermi atoms prepared in two 
different hyperfine states. This allows one to study the 



collective modes across the BCS-BEC crossover region 
in great detail (see, for example, Refs. 0, 0), with ex- 
tremely large values of \a\. One also expects that even 
at the unitarity limit (\a\ — > oo) of trapped Fermi su- 
perfluids, one is dealing with a region where the Landau 
two-fluid equations are now the appropriate description 
of the dynamics at finite temperatures. Several recent 
papers 

nam 

have considered first and second sound 
in a uniform two-component Fermi gas in the unitarity 
region. 

So far, discussions of collective modes in a trapped 
Fermi gas with strong interactions have been limited to 
either T = 0. where one is dealing with a pure super- 
fluid [3 El El EE E 13 El, or to T > T c , where 
one is dea ling with the hydrodynamics of a normal Fermi 
liquid [l!J,|2(J, Um- The much richer collective mode spec- 
trum at intermediate temperatures described by Landau 
two-fluid hydrodynamics awaits exploration. In particu- 
lar, the BEC side of the BCS-BEC crossover is well de- 
scribed as a strongly interacting Bose gas of very stable 
molecules This should be an ideal system to study 
the analogue of first and second sound in a uniform Bose 
superfluid predicted by the Landau two-fluid equations. 

The Landau two-fluid equations describe all super- 
fluids with a two-component order parameter in terms 
of coupled differential equations for the local hydrody- 
namic variables and the velocity fields of the two fluids. 
The coefficients in these equations are given in terms of 
equilibrium thermodynamic functions. Thus, the only 
differences between the hydrodynamics of various super- 
fluids are the values of these thermodynamic functions. 
These depend on the thermal excitations which are ap- 
propriate for different systems (superfluid helium, BCS 
superconductors, trapped Bose gases, etc.). The two- 
fluid equations for a BCS superfluid have been derived 
by Galasiewicz |23j . 

The usual derivation Q of the Landau two-fluid hy- 
drodynamic equations is based on conservation laws and 
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macroscopic considerations, with key equations describ- 
ing the irrotational superfluid velocity field v s (r, t) being 
postulated. There is a considerable literature on this ap- 
proach, which makes little contact with any microscopic 
theory of the origin of the superfluid component. Since 
the 1960s, however, it has generally been accepted that 
this originates from Bose-Einstein condensation (for a re- 
view, see Ref. |24jL 

In the context of trapped atomic gas superfluids, an im- 
portant development was the work of Zaremba, Nikuni, 
and Griffin (ZNG) in 1999 25], who discussed a sim- 
ple microscopic model of a condensate (described by a 
generalized Gross-Pitaevskii equation) coupled to a ther- 
mal gas of atoms moving in a self-consistent Hartree- 
Fock (HF) field (described by a kinetic Boltzmann equa- 
tion generalized to include the Bose statistical factors). 
The ZNG coupled equations for the condensate (super- 
fluid) and non-condensate (normal fluid) have been used 
to discuss the collisionless (see, for example, Ref. |2f| 
and references therein) as well as the collision-dominated 
regions. In the latter region, ZNG showed that in the 
limit that collisions establish complete local equilibrium 
(loth —* 0), their equations can be reduced precisely to 
the Landau two-fluid equations, with all thermodynamic 
functions being given explicitly for a trapped gas of ther- 
mal atoms with a HF spectrum. This formalism is im- 
mediately applicable to a molecular Bose condensate at 
finite T, on the BEC side of a Feshbach resonance 0- 

For a uniform superfluid, the linearized Landau equa- 
tions (without dissipation due to transport processes) are 
solved in the standard texts II IE El m and Predict 
two sound modes. These are the in-phase (first sound) 
and out-of-phase (second sound) oscillations of the two 
fluids, familiar in studies of bulk superfluid helium. The 
mode frequencies have been calculated for a uniform 
weakly-interacting Bose-condensed gas (see Ref. |2|| and 
Appendix C in Ref. |3(J)- For a trapped superfluid, in 
contrast, it is very difficult to directly solve the Landau 
equations for the normal modes of the gas [5E [jHJ . The 
present paper deals with a new approach for calculating 
the frequencies of the normal modes of a trapped super- 
fluid. 

We develop a variational approach of the kind first 
used by ZNG (for a slightly modified version of the two- 
fluid equations reviewed in Appendix B). Our discussion 
is based on an action given in terms of the velocity and 
density fields of the superfluid and normal fluid, first pro- 
posed by Zilsel in 1950 32]. The stationary point of the 
action corresponds to the solutions of the Landau two- 
fluid differential equations. Considering fluctuations of 
this action about an equilibrium state, we introduce sim- 
ple variational ansatzes for the velocity fields. The end 
result is that the two-fluid hydrodynamic modes can be 
described in terms of two coupled harmonic oscillators. 
The effective masses and spring constants of the two os- 
cillators are given explicitly in terms of spatial integrals 
over various static equilibrium thermodynamic functions 
and their derivatives. The latter can be computed in a 



straightforward manner, once we have a specific model 
for the thermal excitations of the trapped superfluid. 

We illustrate our formalism in Section VI by using it to 
derive first and second sound modes in a uniform super- 
fluid as well as the frequencies of the dipole and breathing 
modes of a trapped superfluid. 

There is an extensive older literature [U HE 03 EE 
HE I37L l38l ] on a variational formulation of the Landau 
two-fluid equations, with little or no reference to an un- 
derlying Bose condensate. Indeed, the kind of variational 
formulation we use was originally inspired by the goal of 
giving a clear derivation of the two-fluid equations for 
superfluid helium in the absence of a more microscopic 
theory based on BEC. In contrast, our interest in such a 
variational formulation lies in its usefulness in finding ex- 
plicit solutions of the linearized two-fluid equations cor- 
responding to collective modes in non-uniform trapped 
superfluid gases. However, in Appendix A, we briefly re- 
view some subtle aspects of the variational formulation 
of the two-fluid equations. 

Our discussion of the two-fluid hydrodynamic equa- 
tions using a variational approach is a natural extension 
of the approach recently used to deal with a pure super- 
fluid in trapped gases at T = pi III Hsf . In Section 
VII, we consider the T = limit of our formalism to 
make contact with this "quantum hydrodynamic" the- 
ory, originally pioneered by Pitaevskii and Stringari |39| 
to extend the mean-field Gross-Pitaevskii theory of col- 
lective modes |40j . A hydrodynamic description means 
that a few local variables are sufficient to describe the 
dynamics. While such a description requires rapid colli- 
sions and local equilibrium to be valid in the case of a 
normal fluid, a hydrodynamic description of a superfluid 
is always correct. Thus it is no suprise that the T — > 
limit of two-fluid hydrodynamics gives the correct quan- 
tum description of the pure superfluid (see discussion on 
p. 170 of Ref. 0). 

The discussion in Sections II-IV is fairly technical. 
Readers primarily interested in the final results can im- 
mediately go to Section V and the examples in Section 
VI. Detailed numerical predictions for the hydrodynamic 
mode frequencies at finite T based on our formal expres- 
sions will be presented in future papers. 

II. LANDAU-KHALATNIKOV TWO-FLUID 
THERMODYNAMICS FOR A TRAPPED GAS 

We note that since two-fluid hydrodynamics describes 
a system in local equilibrium, all thermodynamic quan- 
tities we discuss are functions of position and time. This 
dependence will usually be left implicit. Even in static 
equilibrium in the presence of a trapping potential, most 
thermodynamic quantities will still be position depen- 
dent. 

Our variational principle will make use of a Lagrangian 
density of the form 

C = T-U, (1) 
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where T is the kinetic energy density of the fluid and U 
is the internal energy density. Thus, we need to formu- 
late the thermodynamics in terms of the internal energy 
density U, rather than the total energy E = T + U , as 
is normally done in two-fluid hydrodynamics (see, for ex- 
ample, p. 521 in Ref. |23). Our approach will be to ob- 
tain U from the total energy density used by Landau and 
Khalatnikov (LK) to derive the Landau two-fluid (LTF) 
equations 3i]. For a trapped two- fluid system, we define 
the total energy density E as the sum of the kinetic and 
potential energy densities, 



1 2 1 2 
E = 2 PsWs + 2 P " V ™ + + pC/cxt ' 



(2) 



where p n and p s are the normal fluid and superfluid den- 
sities, v„ and v s are the normal fluid and superfluid ve- 
locities, and p = p s + p n is the total density. Also, 



C4xt(r) = -(uj 2 x x 2 + io 2 y y 2 + lo 2 z z 2 ) 



(3) 



is the harmonic trapping potential divided by the particle 
mass, as follows from our use of the mass density p, rather 
than the number density n. 

Following LK, the total energy density Eq(t, t) as mea- 
sured in a frame of reference moving with the local su- 
perfluid velocity is related to E by 



E = E - ^pv 2 - v s • [p„(v„ - v s )[ 
Using the LK identity for local equilibrium 



(4) 



dE = TdS + pdp + (v„ - v„) • d [p n (v n - v s )] , (5) 

with Eqs. J2J an d ®j we obtain the following identity 
for the internal energy density |4l| : 



dU = (jjl- U cyLt )dp + TdS + -(v„ - v s ) 2 dp n . (6) 

In similar fashion, the LK expression for the pressure P — 
d(EoV)/dV including the effects of a trapping potential 
can be recast in terms of the internal energy density, 
giving 

P=-U- P U cxt + TS + pp+ ip n (v„ - v s ) 2 . (7) 

The two definitions given by Eqs. © and J7J can be 
combined to give 



pdp = dP - SdT - p n (v n - v s ) • d(v n - v s 



(8) 



The three thermodynamic identities given by Eqs. 10, 
, and © define all thermodynamic properties that we 
will require in developing our variational approach. We 
note that Eq. ijfjj implies that 



dU" 
9 P /&.,.. 



M = -5- + U, 



(9) 



III. ZILSEL'S HYDRODYNAMICS 

In both classical and quantum mechanics, one can ob- 
tain dynamical equations of motion by equating to zero 
the variation of some action, usually expressed in the 
form A = f dt(T — U), where T and U are the kinetic 
and potential energies, respectively. To derive the LTF 
equations, we take an action of the form j^l 



.4 



d 4 x 



+ ^P nV n ~ U (p, p n ,S)- pU C xt 



(10) 



Here d 4 x = d 3 r dt, while U is the internal energy defined 
by Eq. 10 and £/ cxt is the harmonic trapping potential 
given by Eq. An action similar to Eq. (|l(Jf) was origi- 
nally used by Zilsel in 1950 to derive the LTF equations. 
In this section, we largely follow Zilsel's approach, with 
some differences noted in Appendix A. 

In taking the variation of the action in Eq. I|10|) . the 
variables v n , v s , p, p n , and S will be treated as indepen- 
dent so that, for instance, 



6U = f&U\ 
Sp ~ \dp) S pn 



(11) 



Two important conservation laws not incorporated 
into the action given by Eq. (|10|) are the conservation 
of mass and entropy, described by: 

^ + V- [{p-p n )v s + Pn v n \ =0, (12) 

— + V-(Sv n ) = 0. (13) 

The variation of the action in Eq. I|10|) must be taken 
subject to the constraints given by Eqs. I|12|) and 113|l . 
An additional constraint not employed by Zilsel describes 
conservation of circulation in the normal fluid and allows 
for the possibility of vorticity. This constraint takes the 
form |33l 13(1 137| 



V • (Sr)v n ) = 0, 



(14) 



where r\ depends on r and t. Even though the inclusion 
of this constraint has no effect on the final form of the 
LTF equations, it eliminates the restrictions on the nor- 
mal fluid velocity field that were present in Zilsel's orig- 
inal work. A more complete discussion of this point is 
given in A ppe ndix A. Following the approach pioneered 
by Eckart |42j, this constraint, along with entropy and 
mass conservation, can be incorporated into the varia- 
tional principle by introducing Lagrange multipliers </>, 
a, and 7 (all dependent on r and t) so that the action 
becomes 



.4 



d 4 3 



:(P- Pn)v 2 s + 



2 



U(p,p n ,S) 



4 



{dp 
— + V • [(p - p„)v s + pnVr, 

\ dS « fa > 



{^ + V-, S „v,,) 



(15) 



It is understood that the time and spatial integration in 
Eq. I|15l) is done between two fixed points with the fluctu- 
ations of the variables appearing in the action vanishing 
at both. 

Setting the variation of the action given by Eq. Ijl5|l 
with respect to p, S, rj, v s , and v„ equal to zero, and 
making use of the thermodynamic identities implied by 
Eq. {HI, we obtain the following relations: 



SA 
6p 



2 V * 



(16) 



SA 
SS 



dot _ _ / ' <9 7 _ 



= o, 



SA 
5r) 



-S 



dt 



V 7 = 0, 



id. 
5v., 



= Ps (v s - V<£) = 0, 



(17) 



(18) 



(19) 



and 



SA 

Sv n 



S ( Va + 77 V 7 ) = 0. (20) 



Taking the variation of the action with respect to p n and 
using Eq. (|19f) . one recovers a thermodynamic identity 
already known from Eq. j(JJ|: 



dU\ 1 2 



(21) 



Additionally, taking the variation of the action with re- 
spect to <j), a, and 7 recovers the continuity equation 
given by Eq. I|12|) . as well as the entropy and circulation 
conservation laws given by Eqs. I|13(l and i|14|l . 

Equations fiT ty -ffll jl can be rearranged to yield the use- 
ful expressions 



V0, 



(v„ - v s ) = Va + 77 V7, 



(22) 



(23) 



7 ^ + r ? - = -T--v„.(v n -v s ), (24) 



^7 ■r. 



0, 



and 



30 
dt 



(25) 



(26) 



Equation (1221 shows that the superfluid velocity is irro- 
tational. Furthermore, comparing Eqs. H13[) and (|14|) . we 
find 



drj 



0. 



(27) 



Another relationship that will be helpful in deriving the 
equation of motion for the normal fluid, found from 
Eqs. (ESI, (EEJ, and {HJ, is 

| V7 -gv, = v„x(vx[f(v„-v,)]). (28) 

Taking the time derivative of Eq. H22fl and the gradi- 
ent of Eq. (|26|l . we find the following equation for the 
superfluid velocity: 

dvs _- ( 1 

dt 



P 



(29) 



Rearranging Eq. 0) to obtain an equation for p, and 
making use of the result 



VU = 



dU 
dp 



Vp 



dU 

dp n 



Vp n 



p,S 



dU 
~dS 



VS, 



(/i - U ext ) Vp + X - (v„ - v s ) 2 V Pn + TVS, (30) 



we obtain the expression 
Vp 



I Q 

-VP--VT 

P P 

VU cxt ~^V 
2 p 



(v„ - v s )' 



(31) 



Thus, Eq. (|29(l for the superfluid velocity can be rewritten 
as 

d_ 

dt 



TT7 + V S • V 



-VT 

P 



1 

P 



VP - VU P 



lpn, 
2 p 



(v„ 



(32) 



To determine the analogous equation for the velocity of 
the normal fluid, we take the time-derivative of Eq. 123|) 
and the gradient of Eq. Then, using Eqs. (T2J, 

lfT3jl . JH}, and as well as the identity v • Vv = 
(l/2)V(v 2 ) —v x (V x v), after some labourious algebra, 
we find 



dt + ^- V 



= -^VT 



1 



-VP - vu c 



p 
r 



Pn p 

2 P 



(v„ - v s y 



V.): 



(33) 
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where the "source function" T is defined by 

r^^i + V-( P „v n ). (34) 

From the continuity equation given by Eq. i|12[) . Eq. I|34|) 
implies that 

dp s 



dt 



v ■ (p s v s ) = -r. 



(35) 



Combining Eqs. (|32[l - l|35|) to obtain an equation of mo- 
tion for the current, j = p s w s + /9 rl v„, the source terms 
of the two components cancel and we find 

d] 

tt: = -VP - pVC/ cxt - /) s v s ■ Vv s - p„v n ■ Vv„ 
at 

-v s V ■ (p s v s ) - v„V • (/9„v„). (36) 

More familiarly, in component form, this equation can be 
written as Q 

= ~~Q^~ [""ij + PsVsiVsj + PnVniVnj] 



dXi 



(37) 



where the index j is summed over. 

Taken together, Eqs. (JUJl, (Q3J), US), and §7$ con- 
stitute the Landau two-fluid equations in the non- 
dissipative limit, generalized to include a static external 
potential. 

In Ref. 32], Zilscl claims that the source term appear- 
ing in Eq. (133(1 does not appear in the LTF equations 
of Landau |2j. However, Landau only derived equations 
of motion for the superfluid velocity, as well as the lo- 
cal current density. These are identical to Eqs. I|29(l and 
(I37jl . apart from the external potential we are including. 
If one works backwards from these two equations to de- 
rive an equation for the normal fluid velocity, the source 
term in Eq. ifffg)) does appear. Thus, the LTF equations 
do include a term that contains T. We also note that the 
appearance of T in the equation of motion for the nor- 
mal fluid velocity is consistent with the results derived 
by ZNG from a microscopic model (for further discus- 
sion, see Appendix B). 

Even though our variational method has been devel- 
oped so that we do not have to solve the linearized LTF 
equations, to facilitate comparison with that approach 
we derive those equations here. For the case where 
v s o = v„o = 0, the linearized LTF equations for the 
velocity fields are 



^ = ^V6T-±V5P- 

dt po po 

= -VSn 



VPo, 



(38) 



and 



dt 



PnO PO 

-V5p 



PnO 



Po 
VST. 



(39) 



Here we have made use of the fact that T is zero in equi- 
librium [2{|. The last term in Eq. (|33|l is second order 
in the fluctuations and thus vanishes in a linearized the- 
ory. We have also made use of the fact that in equilib- 
rium VTo = 0. Furthermore, in equilibrium we also have 
V/io = so that for v„o = v s o — 0, Eq. (|31|l yields the 
well-known relation 



VC/ cxt = - — VP . 
Po 



(40) 



Using this, the linearized hydrodynamic velocity equa- 
tions reduce to 

^£ = ^L VS T -VSP ^±5p, (41) 
ot po po po 



<9v« 
dt 



PnO PO PO 

Combining these, we obtain 



Po 



SpVU c: 



k -Sp. (42) 



(43) 



IV. 



ACTION FOR LINEARIZED LANDAU 
TWO-FLUID HYDRODYNAMICS 



As was shown in the previous section, variation of the 
action given by Eq. (|15|l with respect to p, S, ?/, v s , v„, <f>, 
a, and 7 leads to the non-dissipative LTF hydrodynamic 
equations 32] . To determine the low-energy collective 
modes of this system given by the solutions of the lin- 
earized hydrodynamic equations, we Taylor-expand the 
action about the equilibrium values of these variables. In 
principle, one could then take the variation of the result- 
ing action with respect to the fluctuations of its variables, 
i.e., 5p, SS, (57/, (5v s , (5v„, <$</>, 5a, and ^7, to obtain the 
linearized hydrodynamic equations. These could then be 
solved to find the collective mode frequencies. In prac- 
tice, rather than dealing with the Lagrange multipliers 0, 
a, and 7 in Eq. I|15|) . it is more convenient to introduce 
displacement fields for the two velocities. This allows 
one to incorporate the conservation laws directly into ex- 
pressions for 8p and SS, thereby eliminating the need 
for Lagrange multipliers. We will employ a simplified 
Rayleigh-Ritz method in conjunction with our variational 
approach to obtain estimates of the collective mode fre- 
quencies by using physically reasonable ansatzes for the 
displacement fields of the superfluid and normal fluid. 

Recalling that the conservation of circulation con- 
straint does not affect the form of the equations resulting 
from the variational principle, we will omit this term in 
the action for sake of clarity. However, if we were inter- 
ested in collective oscillations of the normal fluid with 
non-zero circulation, we would need to incorporate this 
term and allow for fluctuations of 77 and 7. 

Considering fluctuations of the action about equilib- 
rium we set p = po + 5p, p n = p n0 + 5p n , S = S + SS, 
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4> = 0o + 8<f>, v„ = v„ + <5v„, v s = v s0 + 8v s , and 
expand the action given by Eq. (|15|l up to quadratic 
order in these fluctuations. As in the above, we as- 
sume that the two fluid components are stationary in 
equilibrium, so that 5v n = v„ and <5v s = v s . Writing 
A = + A^ + A^ H , the contribution to the ac- 
tion from terms which are quadratic in the fluctuations 



A {2) = 



d 4 x 



2^0V S + -PnOV„ 



1 fdfX 



f) s j ssp "KS ,,, 

+0o I V • (5p - Sp n )v s + 5p n v„ | 
+a {V-(<5Sv„)} 



2 \ d Pss, Pn 
(SSf 



(44) 



where we have made use of Eq. (jSJ to rewrite the co- 
efficients in terms of thermodynamic derivatives. Note 
that there is an alternative expression for the coefficient 
multiplying SSSp, given by the Maxwell relation 



(dT/dp) s = (dp/8S) p . 



(45) 



Contributions to the fluctuations of the internal energy 
arising from fluctuations of p n have vanished since, as can 
be seen from Eq. 1)21(1 . the equilibrium value (dU/dp n )s,p 
is zero when v„o = v„o. We do not consider the contri- 
butions which are linear in fluctuations since we know 
from the stationary action principle that A^> = 0. 

Since v s0 = 0, Eq. ij^ gives us V0o = 0. Also, with 
v n o = v s0 = 0, from Eq. lll'-'il) we find Vo?o = since 
r\ is an independent and arbitrary variable. With these 
relations, having integrated by parts the terms in Eq. I)44|) 
which contain </>o and ao, the quadratic action simplifies 
to 



! . ) 1 „ „2 , 1 „ ...! i ( <_^_\ i ,\ , ,\'2 



d X i- Ps ov s + -p n0 v n - (*p) 



-S<f> 



d PJs, Pn 2 \ ds 

98p _ , , 
— + V • {p s0 v s + p„ v„) 
ot 



(ssy 



+Sa 



dSS 
dt 



V • (S Q v n ) 



(46) 



In the present form, the variational principle has been 
reduced to equating to zero the variation of the action 



given by the first two lines in Eq. (|46|l . subject to the 
constraints given by the linearized continuity and entropy 
conservation equations, 



98p _ , > 

+ V • (PsOVg + PnQVn) = 



and 



dSS 
~~dt 



V • (5 v„) = 0. 



(47) 



(48) 



For our purposes, it is more convenient to incorporate 
these two constraints into expressions for dp and SS than 
to use the Lagrange multipliers 8<p and 5a. Introducing 
the displacement fields 

v.(r,t) S Vn M) = ^M, (49) 



dt 



dt 



the linearized continuity and entropy conservation equa- 
tions can be expressed in terms of these fields as 



8p(r, t) = -V ■ [p s0 (r)u s (r, t) + p„ (r)u„(r, t)] (50) 



and 



SS{r,t) = -V-[S (r)Un(r,t)]. 



(51) 



Removing the terms that contain the Lagrange multipli- 
ers from Eq. l(4^l) and substituting these expressions for 
dp and 8S, we find our central result for the action that 
is second order in the fluctuations u s and u„, 



A (2 



\ ( 7T ) [V ■ (p s0 u s + p„ u„)] 2 
2 \ d PJs, Pn 

fV-(Sou 



d PJs, P „ 



(SoUn)] [V • {psQU s + p n0 U n )} 



(52) 



To make our expressions as compact as possible, in 
what follows we will drop the p n term constraining the 
thermodynamic derivatives in Eq. (|52l) . However, it 
should be remembered that in evaluating the various lo- 
cal thermodynamic derivatives, p n is always fixed. 



V. HYDRODYNAMIC MODES DESCRIBING 
FLUCTUATIONS 

The fluctuating quantities 8<fi, 8a, 8S, 5p, v s , and v„ 
appearing in Eq. (|4*BT , have been replaced in Eq. lj5^1) by 
the two displacement fields u„ and u s . Thus, the varia- 
tional principle has been reduced to taking the variation 
of the quadratic action given by Eq. (|52|l with respect to 
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only u„ and u s . The linearized LTF hydro-dynamic equa- 
tions and hence, low-energy collective modes of the sys- 
tem, are completely determined by the variational equa- 
tions 



SAW 
Su s (r,t) 



= 0, 



SAW 
5u„(r,t) 



= 0. 



(53) 



Using Eq. 1(52(1 as well as the identities given by Eqs. 150(1 
and J5TJ, Eq. gives 



lis = -V 



:*(*). 



ss 



dp 
dS 



(54) 



and 



-V 

PnO 



Sp 



SP ( dp , s 



SS 



dT 
~dS 



(55) 



where we have used the Maxwell relation given in 
Eq. gSJ). Taking £F = 8S(dF/dS) p +Sp(dF/dp) s , where 
F denotes either one of p or T, with the displacement 
fields defined by Eqs. (|49|l . it is apparent that Eqs. 1(54(1 
and (|55|l are equivalent to the linearized two-fluid hydro- 
dynamic equations given by Eqs. 1(38(1 and 1(39(1 . This 
verifies that the variational principle given by Eqs. i(52|) 
and 1(53(1 is the correct one. 

Solutions of Eq. 1(53(1 corresponding to collective modes 
will have harmonic time-dependence of the form 

u s (r,t) = u s (r)cos(wi), u„(r,<) = u„(r) cos(wt). (56) 

Because these are exact solutions of the variational equa- 
tions, we can insert these expressions into the action di- 
rectly. Substituting Eq. 1(56(1 into Eq. 1(521) and perform- 
ing the time-integration, we obtain, apart from an irrel- 
evant constant factor, the following Lagrangian: 



where 



L (2) = K [u s , u„]w 2 - U[u a , u r , 



K[u a ,Un] = - J d 3 r {p s0 u 2 s + p n0 ul} 



(57) 



(58) 



and 



U[u s ,u 




— [V • (psOUs + p„ U n )] Z 



[V ■ (5 U„)] [V ■ {p s0 U s + PnQUr, 



v ■ (s u n )Y 



(59) 



Since the fields u s (r) and u„(r) are no longer time- 
dependent, it suffices to consider variations of the La- 
grangian given by Eq. l(5Tjl , and the variational equations 
now become 



SLW 
5u s (r) 



Su n (r) 



(60) 



Solving Eq. 1(60(1 is still equivalent to solving the lin- 
earized hydrodynamic equations, with the only simpli- 
fication being that we have assumed harmonic time- 
dependence for the solutions. Fortunately, one can easily 
obtain approximate expressions for the collective mode 
frequencies within a variational approach. Following 
Ref. [25|], we use a simplified Rayleigh-Ritz method and 
make an ansatz for each Cartesian component of the dis- 
placement fields of the form 



i/i0), 



iffi(r), (61) 



where the constant coefficients a S i and a n i are variational 
parameters. Substituting this ansatz into Eq. ((57(1 . and 
equating to zero the variation of the resulting expression 
with respect to these parameters, we have the six varia- 
tional equations 



SIP 
Sa„, 



= 0, 



SJP 
Sa n . 



= 0. 



(62) 



In practice, the symmetry of the problem usually allows 
us to reduce the number of equations. From these equa- 
tions, one obtains a rigorous upper bound for the collec- 
tive mode frequencies u 

A more precise Rayleigh-Ritz method would require us 
to expand each component of the displacement fields in 
terms of n (> 1) elements of some complete set of func- 
tions. As n is increased, solving the resulting 6n vari- 
ational equations, one iteratively generates successively 
better approximations for the collective mode frequen- 
cies. For harmonically-confined gases, however, there ex- 
ist simple trial functions for /j(r) and <7i(r) which are 
sufficiently close to the exact solutions that excellent re- 
sults are obtained by considering only a single expansion 
term as in Eq. (|FTT1) 25] . Our choices of ansatzes for the 
displacement fields at finite temperatures will be guided 
by the known exact solutions at T = ^(j an d T > T c 
(for bosons, see Ref. 0]; for fermions, see Ref. [l^) for 
modes of experimental interest, such as the dipole and 
breathing modes. 

Using the ansatz given by Eq. ((61|) . the Lagrangian 
given by Eqs. I(57(l - I(59|) describes the dynamics of a pair 
of coupled harmonic oscillators, with a S i and a n i repre- 
senting the displacements of the two oscillators from equi- 
librium. The effective spring constants are determined 
by the equilibrium thermodynamic quantities of the sys- 
tem. This is a useful picture to have when envisioning 
the low-energy dynamics of the two fluids. It immedi- 
ately implies, for instance, the existence of in-phase as 
well as out-of-phase oscillation modes of the two fluids. 
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VI. EXAMPLES OF COLLECTIVE MODES 



Uniform Gas 



the following identities which are found from the trans- 
formation properties of thermodynamic derivatives (see, 
for example, chap. 16 in Ref. |4fij): 



As a simple example of our formalism, we consider first 
a uniform gas. In this case, p s o,PnO, and the equilibrium 
thermodynamic derivatives appearing in Eqs. (|58|l and 
(|59|l are all independent of position. Anticipating sound 
modes with a dispersion relation of the form u> = uk, 
where u is the sound speed, we use the plane- wave ansatz 



dp 



dT 



dpj s 2(7 V dp 



dT\ 



a 2 (dT 



dp ) 



1 fdP 



P \ d(T J P P\dp 



dT\ a /«9T N 
dp) a P\d<r /p 



(68) 



(69) 



u s (r) = x J\fa s cos(ujx/u), 



u„(r) = x Na n cos(ujx/u), 



(63) 



where a s and a n are the variational parameters. The 
normalization constant M is chosen so that J d 3 r u 2 = 
a 2 .. Inserting Eq. into Eqs. (SHI and J2I|), we find 



K[a s ,a n ] = ^Psoaj + TjPnoal 



(64) 



and 



U[a s ,a n ] = ^ y 



W ( 




dp 

SqPso 



dT 
dp 



(65) 



The trial functions given by Eq. I|63[l are, in fact, exact. 
Thus, the variational equations given by Eq. which 
in the present case take the form 



0. 



SIP 
8a n 



0, 



(66) 



will give the exact eigenvalues u> or, equivalently, the ex- 
act sound speeds u. Using the above expressions for K 
and U, Eqs. lf5T|) and l(6^|) give the following quadratic 
equation for u 2 : 



dpJs \dp /s 



(So) 2 

PnO 



dT 

~b~s 



Pso , Q , 2 d{p,T) 
WoJ ~tt7 — Frr 

p n0 0(p, S) 



= 0. 



(67) 



Here d(p,T)/d(p,S) denotes the Jacobian of the trans- 
formation between p, T and p, S. 

In the usual textbook discussions, one works with 
the entropy density a — S/p, rather than the entropy 
5 To express Eq. lj6T|) in terms of a, we use 



d(p,T) 1 fdP 
d(p,S) p j 



2 ' dp J rr \ da J , 



dT\ 



and 



dT 



i (dT 



p P\dv /p 



(70) 



(71) 



With these relations, Eq. (|67|l for the sound speed can 
be rewritten as 



dp ) a Pno ° \da 



Pso 2 ( dP 
<X 



PnO 



dp 



dT 



0. 



(72) 



This is the standard equation giving the first and second 
sound speeds in a uniform superfluid @, 0, IH El • 



B. Dipole Mode 

We next consider a dipole mode in a harmonically con- 
fined gas. This mode is characterized by displacements 
of the centre-of-masses of the two fluids along one of the 
axes of the harmonic trap, Xi. In this case, we use the 
following ansatz for the displacement fields: 



(73) 



where a s and a n describe the displacements of the centre- 
of-masses of the two fluids from the trap centre. Substi- 
tuting this ansatz into Eqs. (|58|l and (|59|) . we find 



K[a S) a n ] = ^M s a 2 s + ]^M n a 2 n 



(74) 



and 



U[a s , On] = 7^k s a 2 s + ^k n a 2 + ^k sn (a s - a n f , (75) 

where M s and M n are the the masses of the superfluid 
and normal fluids, respectively, given by 



M s = / d 3 r p s0 , M n = / d 3 r p n0 



(76) 
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The spring constants k s ,k n , and k sn are given by 

dpso 



ks — 



d 3 r 



dp\ dpa ( ( dp\ dSo 
dp J s dxi 



dS J p dx 



d A r 



dji\ dpo 
dp J s dxi 



dp 

dS 




dT y 


) dp ° + 


fdT^ 


\ dS ~ 


dS 


dp j 


1 s dxi 


[dSy 


' P 9xi _ 


dxi 



(78) 



and 



, , dfi\ dp n0 dp s0 , Z' 9T°\ dS dp s0 



dp J s dxi dxi \dp J s dxi dxi f 

(79) 

The variational equation, equivalent to the expression for 
the homogeneous gas given by Eq. if^ with K and U 
given by Eqs. (|71|) and J7SJ), gives the following secular 
equation for the dipole frequencies: 



M s lo 2 - k s 

ksn 



k 

M n aj 2 — k n 



= 0. 



(80) 



Applying some thermodynamic identities, the spring 
constants k s and k n can be simplified considerably. 
Specifically, with v n0 = v s0 = 0, Eq. © can be used to 
find useful expressions for the gradient of various equilib- 
rium thermodynamic quantities. One finds, for instance, 
from the definition T = (dU/dS) p , 



VT 



/ d 2 u \ 

\dSdp) 



Vpo 



Similarly, using Eq. ©, 



dp 
dp~ 



Vp 



dp 
dS 



[ds^) p VSo 



VS + VC/ c: 



(81) 



(82) 



which, it can be shown, is equivalent to Eq. eval- 
uated at equilibrium, when v„o = ~v s q = 0. For a har- 
monic trapping potential given by Eq. noting that 
V/io = 0, Eq. I|82|) provides us with the useful expres- 
sion 



dp \ dpo ( dp\ dSo 



dp J s dxi 



dS j p dxi 



_d_ 

dxi 



■U c: 



-Wi x. 



(83) 



where uji is the trap frequency along the x^-axis. Also, 
since VT = 0, Eq. JS1J gives us 



dT 
dp 



dpo 

dxi 



dT\ dSv =Q 



dS J „ dxj 
p 



(84) 



Substituting Eqs. (jHSJ and into Eqs. and JZHl 
and integrating by parts, the expressions for k s and k n 
readily simplify, reducing to k s = w 2 M s and k n = uifM n . 
With these values, Eq. l{8"01l becomes 

[M s M n (uj 2 - to 2 ) - k sn (M s + M n )] (uj 2 - uj 2 ) = 0. 



(85) 



(86) 



This gives an in-phase dipole mode with frequency 

UJ = LU l 

and an out-of-phase dipole mode with frequency 

^ =^ + ^m n ~ ksn - (87) 

The result given by Eq. I|86fl is the expected Kohn 
mode in a harmonic trap and corresponds to the solu- 
tion a s = a n . This mode is a rigid oscillation of the 
superfluid and normal fluid static distributions and as 
a result, the interactions have no effect on the collec- 
tive mode frequency. The frequency of the second mode 
given by Eq. (|87|l does depend on interactions as these 
determine the thermodynamic functions appearing in the 
k sn spring constant given by Eq. (|79|l . This mode corre- 
sponds to the solution AI s a s + M n a„ = 0. This is pre- 
cisely the analogue in a trapped system of second sound, 
since the displacements of the superfluid and normal fluid 
have opposite signs, producing an out-of-phase oscillation 
of the two components. 

The expression for the frequency of the out-of-phase 
dipole mode given by Eq. (|87|l agrees with that found in 
Refs. |2^| and 0], but with a slightly different value of 
k sn in Eq. I|79|l as a result of working with a modified form 
of two-fluid hydrodynamics (see Appendix B for further 
discussion) . 



C. Breathing Modes 

For the breathing modes, we consider an ansatz of the 
form 



u s = (a sl x, a s2 y, a s3 z), u„ = (a nl x, a n2 y, a n3 z) 
Substituting Eq. JSHJ into Eq. ifS"?). we find 

K[a s , a n ] = \ [ M « a li + M ™°™] ' ( 89 ) 



and 



(90) 
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where the weighted masses Mj are defined by 

M si = / d 3 r PsQ xl M m = / d 3 rp nQ x^ (91) 



and the constants k 



related to spring constants) are 



and k sn< ij (which can be 



dp \ d(p s0 Xi) d(p s oXj) 



dxi 



dxj 



(92) 



dp,\ d(pnQXj) d(pnQXj) 




dxi dxj 



? '9T\ d(pnoXj) d(SoXj) 
dp J s dx i dxj 



dT^ d(S Xi) djSpX j, 
dS 



dxi 



dxj 



(93) 



and 



) fdp,\ d(p s oXi) d(pnOXj) 



dp J s dxi dxj 

dT\ d(p s0 Xi) djSpXj) \ 
dp J s dxi dxj | 



(94) 



Since we generally have six variational parameters (one 
for each Cartesian component of the two displace- 
ment fields), the collective mode frequencies arc found 
from the six coupled algebraic equations, SL^/Sa S i — 
0, 6LW/6a ni = 0. We find 



M si u> a S i = 
2 [( fcs >y 



M ni io a 
2 ^ 

3 



(95) 



(96) 



In experiments performed on trapped superfluid gases, 
one usually has an axisymmetric trap such that, for in- 
stance, uj\ == LO2 = WX; ^3 = w«i and the modes of 
interest are the axial and longitudinal breathing modes 
(see, for example, Refs. |8j and [9j which deal with Fermi 
gases close to unitarity). These correspond to o s i = a S 2 
and a n i = a n2 . 



VII. HYDRODYNAMIC THEORY AT T=0 

Since the variational principle given by Eqs. JS7J, (|61|l . 
and l|t)2|l involves some fairly complex expressions, it is 



useful to consider the case of zero temperature where the 
Lagrangian simplifies tremendously. At T = 0, the Lan- 
dau two-fluid equations reduce to a single hydrodynamic 
differential equation that was first used by Pitaevskii and 
Stringari (3^ to derive the corrections to the mean-field 
Gross-Pitaevskii results jijj for the collective mode fre- 
quencies in a dilute Bose gas. 

At T = 0, only the superfluid component persists and 
we have S = 0, p n o = 0, and p s o — po- Eqs. if5Tjl - if59"|) 
then reduce to 



2, ,2 



(97) 



From SL^/Su s = 0, one obtains 

'dp, 



U! 2 U, 



dp 



V • (poUs) 



(98) 



Using the T = linearized continuity equation, 5p 
— V • (pou s ), Eq. l|9"5|) can be rewritten as 



dp 
dp 



Sp 



(99) 



Multiplying both sides of this expression by po and tak- 
ing the divergence, we obtain a closed equation for the 
density fluctuations, 



u> 8p - 



dp 
dp 



bp 



(100) 



that is the basis of the T — hydrodynamic theory devel- 
oped by Pitaevskii and Stringari [3!j • Eq. I|100|) describes 
the low-energy collective modes of both atomic Bose and 
two-component Fermi superfluid gases at T = 0. The 
only difference between the two quantum gases lies in 
the equation of state p(p). 

It is evident that our variational approach will give ex- 
pressions for the collective mode frequencies at T = 
that agree with results derived from solving Eq. \ 100(1 
directly (see, for example, Ref. |3^| for Bose gases, and 
Refs. m m m for Fermi gases close to unitarity) . How- 
ever, it is still useful to show cxplictly how our formalism 
reproduces the T — results obtained in the recent liter- 
ature. As a specific application, we consider the breath- 
ing modes which were discussed at finite T in the pre- 
vious section. The breathing modes of trapped Fermi 
gases close to unitarit y h ave been studied extensively at 

:r=ommQHEimm, where a simpic ansatz 

for the density dependence of the chemical potential p{p) 
allows one to obtain simple expressions for the frequen- 
cies of these modes. We show that in the limit of zero 
temperature, our general results for the breathing mode 
frequencies at finite T, given by Eqs. (|95|l and H96fl . re- 
produce these well-known expressions. 
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Since the normal fluid vanishes at zero temperature, we 
only have equations for the superfluid component, given 
by Eq. (|95(l . which reduce to 



M si uj 2 a si =y^k s 



(101) 



at T = Q, where k n and k sn are zero. The expression for 
k s ,ij, given by Eq. ll§2*|) . becomes 

fcs „ = /d 3 r (dt) d( P0Xi ) d( P0Xj ) ^ 



dp J dxi dxj 



since p s o — po at T = 0. 

To obtain a simple analytic expression for the collec- 
tive mode frequencies, the following identity, found by 
integration by parts, will prove useful: 



dp\ dpo 



dr [ dp-) P ° Xj d^ 



dfi 



d [dp 



-2j drp i{t) +x ^{fp)j- (103) 

Expanding out the derivatives in Eq. H102J) and using 
Eq. I|103|) . Eq. (|1U2[I can be rewritten as 



k,, = I d 3 r 9P0 9P0 



dp ) dxi dx.j 



Xi X j 



2 \ dxj dxi ) \ dp 



At T = 0, Eq. (|H31) reduces to 



dp\ dp 



004) 



(105) 



dp ) dxi 



Using this result in Eq. i|l(J4|) and integrating by parts, 
we find 



k s , t j = J d A r po j(2% + lMxf 

po d ( dp\ po d ( dp" 



2 3 dxj \dpj 2~ l dxi \dp / 



(106) 



Assuming a so-called yolytroyic equation of state (see, 
for example, Refs. EH El 13 El El ) , 



we find 



H[p] tx p 7 , 



a dxi\dpJ \dp J dxi 



(107) 



-(7 - 1K 2 ^ 



(108) 



Substituting this relation into Eq. (|106fl . an d making 
use of the definition of the weighted mass M S i given by 



Eq. (|91|l , we obtain the following expression for the ratio 
of the spring constant and the weighted mass: 



^3 



2uj 2 5 lj 



(7 + 1) 



+ Xji 



(7-1). 



where 



Xji 



J d 3 r p (r)x 



Jd 3 r p (r)x 2 



(109) 



(110) 



Using the Thomas-Fermi result po = p[po( r )] + U cx t(r), 
the equilibrium density po can be expressed as a function 
ofC/ C xt(r) = {1/2){uj\x\-\-u>2x\-\-lo1x\). Making a change 
of variables, x' = uJiXi, \ji becomes 



w? / d 3 r' p [xf + x' 2 2 + x'i)x' 



J d 3 r' p [xf + xf + xf]x: 

, ,2 
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With this result, Eq. (|109|) simplifies to 

ks,ij o 2 r i 2 

Combining this identity with Eq. I|101|l . we find 



(HI) 



(112) 



u 2 a x 



2a; 2 a, 



-7K 



(113) 



which is equivalent to Eq. (3) in Ref. |l8j when a poly- 
tropic equation of state is used. 

For an axisymmetric trap, u>i = uj-2 = cjxj ^3 = w 2 , 
and the axial and longitudinal breathing modes are char- 
acterized by solutions of the form a\ = = a. In this 
case, Eq. (|113fl gives the normal mode frequencies found 
by previous authors pi. HI UlEsT] . 



2(7 + lK + ( 7 + 2K 2 



v/[2( 7 + l)u£ - (7 + 2K 2 ] 2 + S^ui 



(114) 



VIII. CONCLUSIONS 



More than sixty years after Tisza Q] and Landau 
proposed two-fluid hydrodynamics, their theory turns 
out to be very relevant in the study of the dynamics 
of trapped Fermi superfluid gases with strong interac- 
tions, close to unitarity. From linearized two-fluid hydro- 
dynamics, one can obtain collective mode frequencies in 
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terms of equilibrium thermodynamic functions and ther- 
modynamic derivatives. This is well known for uniform 
superfluids and still occurs in the case of trapped gases, 
apart from the fact that the equilibrium thermodynamic 
quantities are now dependent on position. This means 
that the two-fluid predictions are general, with a micro- 
scopic theory only being needed for the elementary ex- 
citations that determine the equilibrium thermodynamic 
functions. 

Variational approaches have already been used exten- 
sively to obtain hydrodynamic collective mode frequen- 
cies of superfluid Fermi gases at T = close to unitar- 
ity 0, 0, ^3 1 as we ^ as boson-fermion mixtures • 
In this paper, we have extended this approach to finite 
temperatures to deal with the coupled hydrodynamics of 
a superfluid and normal fluid in a trapped gas. 

We have presented a variational method to determine 
the collective mode frequencies of a trapped superfluid, 
based on an action of the type first proposed by Zilsel 
in 1950 to derive Landau two-fluid hydrodynamics. By 
making an ansatz for the superfluid and normal fluid 
velocity fields, in this formulation, the problem of de- 
terming the collective mode frequencies is reduced to 
solving a set of algebraic equations. This obviates the 
need to solve the Landau two-fluid differential equations 
numerically, a difficult task in a nonuniform trapped gas. 

We illustrated our formalism by considering a uniform 
gas, as well as the dipole and breathing modes for a har- 
monically confined gas. To make contact with recent 
discussions in the literature on trapped superfluid Fermi 
gases, we also used our finite T variational approach to 
discuss the collective modes in the limiting case of a pure 
superfluid at T = EHU. 

In future work, we will apply the expressions derived 
in this paper and present detailed numerical predictions 
for the frequencies of various hydrodynamic modes in 
trapped superfluid gases at finite temperatures. In par- 
ticular, we will discuss results in a Fermi superfluid for a 
strongly interacting molecular BEC at finite T (the BEC 
side of a Feshbach resonance) . 

In this paper, we have only discussed how to determine 
the normal modes of the non-dissipative Landau two-fluid 
equations for a trapped (non-uniform) superfluid. The 
damping of these modes due to various hydrodynamic 
relaxation processes can be calculated using the general 
expressions derived by Nikuni and Griffin |50|. 



Acknowledgments 

We would like to thank Duncan O'Dell for stimulating 
discussions about his work on hydrodynamic modes in 
the unitarity region of Fermi superfluids. We also thank 
Eugene Zaremba for useful comments related to the vari- 
ational approach developed in Ref. |25|. Our research 
was supported by a grant from NSERC of Canada. 



APPENDIX A: COMPARISON WITH ZILSEL'S 
VARIATIONAL PRINCIPLE 

As discussed in the text, our major interest in this pa- 
per is to find an efficient way of calculating the collective 
mode frequencies of a trapped superfluid in the region 
where Landau's two-fluid hydrodynamics are valid. Fol- 
lowing Zilsel [33], we did this in Sections II-IV by refor- 
mulating the two-fluid equations in terms of a variational 
theory involving the action. There is an extensive older 
literature on such a variational formulation in the con- 
text of superfluid helium. In this Appendix, we review in 
more detail some of the subtle issues which arose in con- 
nection with Zilsel's seminal paper and their resolution. 

One criticism concerned Zilsel's use of a; = p n /p as 
an independent variable which, it was argued, conflicted 
with the Landau formalism 34]. Zilsel was ultimately 
vindicated by Jackson who demonstrated that Zilsel's ap- 
proach was compl etely equivalent to those proposed by 
some of his critics [3(j . 

In our approach to the two-fluid thermodynamics, us- 
ing the internal energy density given by Eq. ©, we have 
deviated slightly from the variational method devised by 
Zilsel to derive the LTF equations. He considered the 
specific internal energy, which we denote by Uz to avoid 
confusion. The two energies are related by U = pUz- 
Furthermore, while we have taken the internal energy 
density to be a function of the entropy S, the total den- 
sity (p = p n + p s ), and the normal fluid density p n , Zilsel 
took the specific internal energy to be a function of the 
entropy density a — S/p, p, and the ratio x = p n /p- 
However, one can show that the thermodynamics given 
in Section II, obtained from the Landau-Khalatnikov en- 
ergy density, is identical to Zilsel's and our choice of in- 
dependent variables is completely consistent with Zilsel's 
choice. To compare our thermodynamic identities with 
Zilsel's, using Eqs. © and J7J, with U = pUz and hence, 
dU = pdUz + Uzdp, we find 

P 1 

dU z = Tda+— dp+ -(v n -v s ) 2 dx, (Al) 
p z I 

where a — S/p and x = p n /p. This relation gives pre- 
cisely the thermodynamic relations used by Zilsel in his 
derivation of the LTF equations. As well, Eq. I|A1(I was 
derived by Jackson using thermodynamic arguments |3fi| | . 
Zilsel chose to use ct, p, and x as independent variables 
as these differentials appear in Eq. (|A1|) . In contrast, 
with the internal energy density given by Eq. JfJJl, the 
differentials of S, p, and p n are the natural choice for 
independent variables. 

A more valid criticism directed at Zilsel's work con- 
cerned the fact that, while he obtained the full LTF equa- 
tions, his variational principle forced the normal fluid ve- 
locity to be irrotational when S/ p n is a constant. Thus, 
Zilsel's variational principle only describes a subset of all 
possible solutions of the LTF equations. Following an 
argument originally due to Lin, this deficiency can be 
corrected by introducing an additional constraint, given 
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by Eq. i|14|l . into the variational principle that accounts 
for conservation of circulation (or vorticity) in the nor- 
mal fluid HI El El El H| To see this explicitly, 
consider Eq. (|23[) . Without including the circulation con- 
straint, equivalent to setting rj = 0, we recover Zilsel's 
expression for the normal fluid velocity. It is apparent 
from Eq. 1|23|) that for the case when the entropy divided 
by the normal fluid density is uniform, the normal fluid 
velocity is irrotational. By including the circulation con- 
straint, the difference between the normal and superfluid 
velocities multiplied by p n /S assumes the form of a Cleb- 
sch transformation [5 If , with the local functions rj and 7 
acting as the Clebsch potentials. In this case, the normal 
fluid velocity is no longer restricted. 



APPENDIX B: ZNG VS. ZGN 
HYDRODYNAMICS 

We briefly summarize some key features of the deriva- 
tion of two-fluid hydrodynamics for a trapped Bose gas. 
As discussed in detail in Refs. [25| and j2ij, the ZNG 
microscopic coupled equations involve in a crucial way 
a term Ti2 which is related to the collisionless transfer 
of atoms between the condensate and thermal cloud. In 
the linearized ZNG formalism, one can show that in local 
equilibrium, 



sr 12 = 



(Bl) 



This involves the difference in the chemical potentials of 
the condensate (/i c ) and the thermal cloud (/i) 



(B2) 



In static equilibrium, Ti2 vanishes and fio = Mco, assum- 
ing that v c o = v„o = 0. In Eq. i|Bl|l . the relaxation time 
T12 describes collisions between atoms in the condensate 
and thermal gas components. The relaxation time r M in- 
troduced in Ref. [25| which determines the rate at which 
local diffusive equilibrium is achieved (fl = /i c ), is given 
by r M = OffTi2, where an is a dimensionless hydrody- 
namic normalization factor defined in Ref. p5l | involving 
local equilibrium thermodynamic functions. 

ZNG derive the non-dissipative Landau two-fluid equa- 
tions in the limit of wt m « 1. In this limit, one finds 



— > 0, but 8T12 in Eq. (|Bip remains finite and is 
given by 

ST 12 (r,t)=a H lv • [n c (v c - v„)] + ^n c V ■ v„ j . (B3) 



As discussed in Ref. j25j, this source term in the conti- 
nuity equations Eqs. (|3"4"|l and (|331) plays a crucial role in 
establishing the equivalence between the ZNG equations 
(in the limit that lot^ — > 0) and the standard Landau 
two- fluid equations considered in the text of this paper. 
Thus, Ti2 is a concrete example of the source term T 
(= TOL12) which arises in Eqs. (|33|l and l|34|l . in the con- 
text of a dilute Bose gas. 

This is a convenient place to note that in Refs. pol 
I29I |46| , another version of the two-fluid hydrodynamics 
was discussed in which the collisional transfer of atoms 
between the condensate and thermal cloud was ignored. 
This corresponds to setting Ti2 = or, cquivalently, 
the limit ujt^ ^ 1. This "ZGN" two-fluid hydrodynam- 
ics was used to calculate the temperature-dependence of 
the frequencies of first and second sound in a weakly- 
interacting uniform Bose gas [23. In this case, the veloc- 
ities were almost identical (see Fig.l of Ref. j25j) to those 
based on the Landau two-fluid hydrodynamics (wr^ <C 1 
limit) . Expressions for the two sound velocities based on 
the LTF equations are given in Appendix C of Ref. [13, 
showing the small corrections which depend on the factor 
(JH in Eq. (|B3)l . 

The variational theory developed in Ref. [2I] for de- 
termining the normal mode oscillations in a trapped 
gas was worked out originally for the ZGN two-fluid 
hydrodynamics [4fij . The results for monopole, dipole, 
and qua drupole oscillations are summarized in Fig. 1 of 
Ref. [46| and Figs. 4 and 5 of Ref. 0. For the dipole 
modes, the frequencies are identical to Eqs. 186|) and 18711 . 
but with 



-2g / d 



dn c p dh 
dxi dxi ' 



(B4) 



where g — Ai:Ti 2 a/m. It can be shown, however, mak- 
ing use of approximations consistent with those used to 
derive ZGN hydrodynamics, that Eq. I|79|l reduces to 
Eq. (|B4|I . The details of this derivation will be given 
elsewhere. 

The variational approach of Ref. [2j| was formally ex- 
tended to deal with the ZNG hydrodynamics in Ref. [5(J , 
but it was not used to calculate the collective mode fre- 
quencies. Instead, the out-of-phase mode was approxi- 
mated by extending the T = pure condensate mode 
to finite T by using a reduced condensate density. The 
in-phase solution was approximated at finite T by the 
thermal cloud oscillation above T c , but with a reduced 
thermal cloud density. As discussed in Section VII of 
Ref. [25j], these are often good first approximations. 
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